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Abstract 

We study in a rigorous way the XYZ spin model by 
Renormalization Group methods. 

pacs 75.10.Jm 

Renormalization Group, Correlations, XYZ 
chain 

The XYZ model is exactly solvable [1] by the 
transfer matrix formalism as it is equivalent [2] to the 
eight vertex model. The solution is so complicated 
that it is very difficult to compute the correlation 
functions from it (an attempt with some preliminary 
results is in [3]). The correlations are then computed 
in an approximate way by linearizing the bands and 
taking the continuum limit [4], so introducing spuri- 
ous u.v. divergences. One has to introduce an ad hoc 
cut-off, absent in the original model, for applying the 
bosonization methods, and it is not very clear the re- 
lationship of the obtained correlations with the real 
ones. Finally Bethe ansatz and the conformal algebra 
methods cannot be applied to the XY Z chain but 
only to its limiting case given by the XX Z chain [5] 
(also bosonization results are mainly for this case). 

In this letter we apply to the XYZ chain the RG 
methods developed for QFT [6]; we show that, for 
small anisotropy and J3, the correlations can be ex- 
pressed by convergent series. This is almost equiva- 
lent to know the correlations exactly, as one can com- 
pute the first orders having a rigorous bound on the 



remainder i.e. the correlations are known up to a 
small error. The comparison of the so obtained corre- 
lations with experiments or numerics is then without 
ambiguity. No approximation are necessary in our ap- 
proach, which can be extended to a variety of models. 
Bounds on the corrections due to the finite size effect 
are naturally obtained, as our results are uniform in 
the size. We write the XYZ chain as a system of 
interacting fermions, and we compute the fcrmionic 
two point different (imaginary) time correlation func- 
tion. The computation of the other correlations is 
a straightforward (but cumbersome) consequence of 
our results. 

If (S^SlS*) = \{al,al,a% a°, a = 1,2,3 be- 
ing the Pauli matrices and x — 1, 2, N the hamil- 
tonian H is 

N 

[J\S X S X+1 + J2S X S X+1 + JsS x S x+1 + hS x ] + U N 



where the last term is a boundary term and S 3 



N+l 



Sf (the boundary conditions on S x , S x will be speci- 
fied later fixing U^)). The anysotropy is u = 



_ ./1-./2 
.h+.h 



> 



for fixing ideas. 

By a Jordan-Wigner transformation [7] we write, 

if 5± = S x ±iSl, S~ = e""^»=i ^~ and S+ = 

ip+e^^y^ Jp ^ Jpy , where ibt 



ip x are fermionic operators. 



Moreover S x = ip x ip x — |. Fixing Ji + J 2 = 1, H is: 



iV 



M^ZA -gHV^+iVwi- 5) +M^V>a - \)} + u 2 N 



1 



S+e^Si) + 



S%e mM Si ) + c.c], if M 



where Ujj is a boundary term. We choose 
so that [/jy = and the fcrmions verify peri- 
odic boundary conditions (in [7] this choice for the 
XY chain is called "c-cyclic"). This is achieved 
by setting - 
c.c.}+ u[(-S+S+ + -L- c+.i^c- 
Y^,x=i Vkj as [(— 1)^,-7?] = the eigenvectors of 
H are divided in two subspaces on which (— l)^ is 
equal to 1 or —1 and on such subspaces is a 
boundary term. Let be S £l,e2 S(x) = lini7v,/3^oo < 
Ttp^tpi 2 >N,p, £i = ±, T is the time ordering and < 

• >= tre-^H • We P rovc that > if hmjv.^oo S^p*(x) = 
S ei ' £2 (x), x = x,t and setting pj? = cos _1 (/i — J 3 ) ^ 
nn, vq — sinpi? and \x\ = \J x 1 + v^x^, k — k, k : 
For J3 , u suitably small 

S ei < E2 (x) = S e ue2 {x) +eSl ue2 {x) 

where Sq' + (x) and Sq' + (x) are respectively 



given below; if 11 = a(k) = and Z(k) = \k\~' n3 
according with the expected power law decay of the 
XXZ chain. A simple consequence of our analysis 
is that, in the J 3 = h = case, | < S^Sq > - < 
52 >< 5l > I < T§ra-e~ C2U l £ l in agreement with the 

x 1 — \ x \ 

results in the XY chain in [7], for small u and t = 
sin 2 (fx)^(^r4,[l + 0(i)]. 

i?G analysis. The Grand-Canonical Schwinger 
function can be defined by Grassman integrals so we 
can write < Tn"=i >£,0 as 



/{^/^'^"^'^sje-^l n"=i $3. 



(1) 



where k 



dk sin(p F x+kot) 
J dksin(p F x 



2ie l 



ik ~ v sin(|fc| - p F ) 



Z(k) k%+v%sm 2 (\k\-p F )+a(k) 
2ie tkS cr(fc) 



Z(k) fcg + w 2 S in 2 (|fc| - p F ) + <j(k) 2 

where e = max(w, | J3I, u l+m ), p F — p F + 0{J^) and 
Z(k),a(k) are smooth functions such that \Z(k) — 1 < 
e, \a(k)-u\ < e for \\k\-p F \ > mm[(p F /2, (>K-p F )/2] 
and 



Z(±p F ) = z = u -' ni 



a(±p F ) — a = u 



= rr — nA+m 



within =Pi{-h) 2 +0({hY), m = ~(3 2 -h + 0{{h) 2 ), 
/3i,/?2 > 0. Moreover for \x\ > vT 1 , Sq 1 ^ 2 ^) have a 
long distance faster than any power decay, i. e. for any 
M, \S £ ^ £2 (x)\ < if c M ls a cons tant; 

for 1 < \x\ < u^ 1 they have a transient slow decay 

\s E ^ 2 (x)\ < with m = mii + m)- 1 . 

The optimal bound for S , q 1 1 ' E2 (x) should be 
\S E ^ 2 {x)\ < ^^e-^" 2 !*!, Cl ,c 2 > constants, 
and this could be proved by a slight improvement 
of our techniques. We see that there is an anoma- 
lous gap and an anomalous wave function renormal- 
ization. Moreover the oscillation period p F depends 
on J3. The explicit form of Z(k) and cr(k) will be 



^(tv' 212 ^^ — )) ^ n o,ni are integers, 
E(k) = cosfc — cosp F , are Grassman vari- 

ables and we denote by J {Vipe f dk ^h h< -^ } the 
fermionic integration with prop- 

agator / dke^-^h(k), if f dk = Finally 
V(i/j) can be written as V(ip) = J3V + uP+vN where 
V = I nti dh cosfo - ^t^l^/ikr -k 2 + 
2 k 3 - fc 4 ) and P = J dk[e ik tpitp+^ + e~ lk tPl^Z], 
N = f dkibtibZ . Note that v is a counterterm to be 

J k k 

fixed so that the Fermi momentum in the J3 = or 
J3 ^ theory are the same. In fact, as the oscillation 
period is changed by the presence of the J3 term, we 
find technically convenient to fix it to its value in the 
J3 = case by adding the counterterm. According to 
the (formal) Luttinger theorem, this means that we 
add a magnetic field so that the mean magnetization 
in the direction of the magnetic field is the same as 
in the J3 = case. 

We start by integrating the denominator of eq. (jl]) , 
the partition function AT. We perform a decompo- 
sition of the propagator g(k) — fi(k))g(k) + (1 — 
h(k))g{k) where /i(fc) = 1 - x(k -pF,h) - x{k + 
p F , ko) and x is a smooth compact support function 
such that x(fc ±p F ,ko) are non vanishing only in 
two non overlapping regions around p Fl and — p F , 
respectively. We call the two addends respectively 
g >0 (k) and g-°(k) and this allows us to represent ipi 
as sum of two independent Grassmanian variables, 
^,(>o)±^ ^(_<o)± rp-^ c m t e g ra tion on -)/A >0 ) allows us 

to write Af with ip replaced by ip 1 -- ^ and V replaced 



2 



by V° defined as [8]: 

P n ii 

j n^^(^,...,4)n4 : 



oo 



(<h)si 



n=0 



£ ih) 



with h — 0, where 
z = max(|A|, \v\). 



(2) 

W°\ < (7 m z max ( 2 ."- 1 ) if 
Note that, unless u = 0, 
ElLi^ ^ as [i?,A/] ^ 0. We set k = k' + uip Fl 
uj = ±1, sign(o;fc) > 0; moreover we write 1 — 

X {j- h+1 k') and C,; 1 = EL-oo h- This allows us to 
write vi- 0) = E^iiEL-oo^- In other words 
we represent the fermions by two Fermi fields with 
label h, with momenta close 0{^ h ) respectively to p F 
or — pp. We proceed iteratively setting Z = 1: once 
the fields . . . , ij;( h+1 ) have been integrated J\f is 
given by: 



{V^e-J d * ChZh ^ h 
-v"(VZ^/> ( -" ) ) 



(fe')- 1 ^7 ,l) 



} 

(3) 



if V- 



I and V h {ip^) 

is given by terms like eq.(H). If a(fc) = cosj5f(1 — 
cosfc), [Q {h) (k')]z^ = {-iko + uvosmk + uja(k')), 
[Q W (k')]~^_ u = ioJa h (k'), We define a focafea- 
izon operator C extracting the relevant part of the 
effective potential V h : i) If n > 4 then £W,J = 
0; ii) if n = 4 if 5f = <5^4 



then CW£(h 



if 5? = 0^4 „U^4 

WiPf,-.) = 5lW%{u>ip F ,...) iii) if 



n = 2 then if <Sf = 8^ 1 ^ UJ2 ) PFfl 8 ei - E2fi and tf£ = 
~(Lu 1 +ui 2 )p F ,oO~e 1 +e 2 ,o then 

+ Vip F ,k' 2 + u) 2 p F ) = 5 2 l [W 2 h {uj 1 p F ,Lu 2 p F ) + 

uj 1 E(k' + uJxPF)dkW 2 h (w 1 p F , u> 2 p F ) + fc° 

d ko W 2 (cJip F ,uj 2 p F )} + 6 2 [W 2 (cJip F ,uj 2 p F )} 

where E(k' + ujp F ) — voujsink' + a(k') (the symbols 
dt, dk means discrete derivatives and the first deltas 
in S 2 ,S 2 ,8\ are mod. 2ir). A naive power counting 
argument explains why the relevant terms are only 
the quartic or bilinear in the fields; moreover among 
such terms there are still irrelevant ones i. e. the power 
counting can be improved. This is taken into account 
in the definition of C, as the first of the two deltas 



in 6%, 5%, S 2 says that the relevant terms involve only 
fermions at the Fermi surface, i.e. J2i £ i UJ iPF = 
modulo 2-7T, and the second takes into account that the 
marginal terms with J^i £ i 7^ are indeed irrelevant. 
This will discussed below. We can write then the 
relevant part of the effective potential as: 

CV h = j h n h F£ + s h F* + z h F* + a h F* + l h F^ 

where F, h = V f dk' fii> ( --. h)+u <ib { ^ h) -" , and 
Fa = E fdk'^ h)+u d,^ h }~ u and J* is 

given by f\nt=i^i^!^*J S ^Li w * tr A) where 
i = u,a,C, fu = uj, f a = LoE(k' + p F ), f c = -ik . 
Moreover l = J 3 + (9((J 3 ) 2 ), s = u + 0{uJ 3 ), 

«0)^o — O(Ja), n o = v + 0(J 3 ) and we have defined 
l p{<h)±u = ^(<o)± if ^ = 1 an d </,(<o)"± = ^,(<0)=F if 

to = — 1. We write eq.(3) as: 



(4) 



where V h = CV h + (1 - C)V h , CV h = j h n h F^ + 
(a h - z h )F* + l h F* and Z h ^(k') = Z h + C^Z h z h , 
Z fc _i(fc')a- h _i(A / ) = Z h a h (k>) + Z h C- l s h . This 
means that we extract from the effective poten- 
tial the terms leading to a mass and wave function 
renormalization. Now one can perform the integra- 
tion respect to ip^ rescaling the effective potential 
V h (iP) = y^(y^Z^) and CV h = j h v h Fl l + 5 h F* + 

\ h F£ with j h u h = ^n h , S h = ^(a h - z h ), 
\ h = (-^^) 2 lh and v k = i/ h ,S h ,\h- The in- 
tegration of ifj^ has propagator g^^/ix — y) = 

^Jdk'e i Z(*-y)f h {k')g( h - l \k') u ^, if z h _ x = 

Z h _i(0) and f h = Z h _ x \ z ^ {kl) - gi]. After this 
integration M is, up to a constant, of the form of 
eq.(3) with h replaced by h — 1, and we can iter- 
ate. In other words we have to perform a Bogol- 
ubov transformation for each scale, as the "mass" 
(T/j has a non trivial RG flow and it is different for 
any h; at the same time one has to take into ac- 
count the wave function renormalization Zh- Let 
be h* = infh{7 /l > |o7i|}- Note that, if h* is fi- 
nite uniformly in N,f3 so that |e>7i*_i|7~' 1 +1 > 1 
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then \g <h *(x)\ < 1+ ^7 |g|)M ; moreover if ft > ft* 

we have \g% tU (x)\ < ^ 1+ f 7 "^ )M and | 5 £ _ u (x)\ < 

-kz^- 1+ ^\ S \)m ■ The propagator for the integration 
of all the scales ft < ft* obeys to the same bound 
of a single scale propagator for ft > ft*. Moreover 
for ft > ft* the bound for the non diagonal propaga- 
tor has a factor more with respect to the diago- 
nal propagator. Finally g„ tU {x - y) = g'^ L (x - y) + 
C^Jx - y) + C^(x - y), with g^ L (x - y) given by 

/ dk'^-g — ^l jfc i^tttffc' ^ which is just the propagator 
"at scale h" of the Luttinger model, and the other 
two terms verify the bound of g^ u {x;y) with an ex- 
tra factor 7 71 or ^rr- 

Lemmavlssume that ft* is is finite uniformly in 
N, (3 and that for any h > k > h* there exists an 
e such that \v h \ < e and \^-\ < Y aE , \^ L \ < 7^ 
with c a ,Cb postive constants. Then V k is given by a 
convergent series. 

We have to show by the study of the beta function 
that the conditions in the lemma are verified. It is 
possible to choose [8], [9] the counterterm v so that 
\vh\ < e for all > ft > ft*. The Beta function can be 
written, for > ft > ft*. : 

x h -i = x h + G\ h + g\ h + ^ h R\ 

o~h-i = o~h + G 1 / 1 

6 h -i = S h + G\ h + Gf h + 7 fc i$ 

^ = l + Gl> h + Gl> h + 1 h R h z (5) 

where a) G^ h , G]' h and G\ h depend only on 
Ao, So; ...Xh, 5h and are given by series of terms in- 
volving only the Luttinger model part of the propa- 
gator g k L (x — y), so they coincide with the Luttinger 
model Beta function [8], [9]; b) G 1 / 1 , G\ h , G 2 ' h , G 2 z h 
are given by a series of terms involving at least a prop- 
agator C% >u {x - y) or g^- u (x - y) with k > ft; c) 
Ri, i — X,z,5 are given by a series of terms involv- 
ing at least a propagator C\ UJ {x — y), k > ft. By 
a simple computation G 1 / 1 = X\ + G% G\' h = 
\ h o h [-p2 +G% with fa, fa > and G h z ,G h a = 
O(Xh). Moreover G 1 ^' 1 , G\ ,h coincide by definition 
with the Luttinger model Beta function, and it was 
proved in [8], [9] that it is vanishing at any order, i.e. 



G\ h {X,5:...;X,5) = and G\' h (X, 5; . . . ; A, 6) = 0. 
Finally as \G^ h \,\G 2 s h \,\G 2 z ' h \ < Ke 2 \o h \T h \ one 
finds, for ft > ft*, |A^_i— A | < c x l%, \b h -\-^\ < CiZg, 

l pic 2 h < log ( Ihzl. ) < ? /3 lC3 ft 
V CT o / 

-Mlh < Iog(Z fc _i) < -(3 3 c 5 l 2 h 

for suitable positive constants Cj, i.e. as usual in 
models to which the RG is succesfully applied the flow 
is essentially described by the second order truncation 
of the beta function. This shows that it is possible to 
choose J3 so small that the conditions of the above 
lemma are fulfilled. We call 771 = — log Z/,. / log it, 

1 +T] 2 = logCT/j./logW. 

Finally as we said the integrations of the ip( <h ) 
is essentially equivalent to the integration of a single 
scale h > ft*, so it is well defined by the preceding 
arguments. 

It is a standard matter to deduce an expansion for 
the correlations from the effective potential, and so 
deducing the results for S ai ' a2 (x). Moreover we call 
a{k),Z{k) respectively ah, Z h for j h < \k\ < ^ h+1 for 
ft > ft* and o- h *,Z h * for \k\ < j h \ 

Finally we discuss the modifications in the proof 
of the above lemma with respect to the one exist- 
ing in literature for similar models, [8], [9], mainly due 
to the fact that [TV, H] ^ 0. V k can be written as 
a sum over "Feynman graphs" obtained in following 
way. Let us consider n points and enclose them into 
a set of clusters v to which a scale h v is associated; 
an inclusion relation is established between the clus- 
ters, in such a way that the innermost clusters are 
the clusters with highest scale, i.e. if v' is the cluster 
containing v then h v > < h v ; v is the largest clus- 
ter. A set of clusters can be represented as a tree and 
the set of the possible trees is denoted by T n ^. To 
each point contained in a cluster v but not in any 
smaller one we associate one of the elements of £V hv , 
expressed graphically as a vertex with 2 or 4 "half 
lines". The half lines are paired in all the possible 
compatible way and to the paired lines we associate 
a propagator g u v u , if the paired line is enclosed in the 
cluster v but not in any smaller one. The indices of 
the external lines of v are denoted by P v and their 
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number by \P V \. The 1Z = 1 — C operation acts on 
the cluster and its action, as we know, depend on P v . 
We call V k the contribution to V k from r G T n ^, and 
\V*\(L(3-y- kD( - p »o))- 1 = \V k \ is bounded by, if Ci is 
a constant 



eik 



E 



[D(P v )+z(P„)](h v -h v ,)l 



where D(P„) = -2 + \P v \/2, and V 2 is the set of 
the clusters v containing a non diagonal propagator 
9v°-ui i- e - the smallest clusters containing a non di- 
agonal propagator. The factor z(P v ) (which would be 
zero if we had naively taken 1Z = 1) is defined as:l) 
z{P v ) = 1 if \P V \ = 4 and 5% = 1; 2) z(P v ) = 2 if 
\P V \ = 2 and 55 = 1; 3) z{P v ) = 1 if \P V \ = 2 and 
8\ = 1; 4)z(P„) = in the remaining cases. The first 
parenthesis is the power counting of a graph, and the 
second is due to the extra factor in the bound for 
the non diagonal propagators; the only non trivial 
part of the above bound is the use of the Gramm- 
Hadamard inequality to take into account the relative 
signs of the graphs (estimating them by their abso- 
lute value one would get an extra fc! in the bounds 
spoling convergence). If D{P V ) + z(P v ) > the 
sum over r, {P v } can be bounded and \V k \ < £ n ; 
however, if 8%, 8%, 62 = this is not the case. The 
factors 5%, 5%, 5% are products of two deltas. If the 
first delta is vanishing, by the support in momen- 
tum space of g^ u /, it follows that there exists a fixed 
scale h, indipendent on L, f3, k, n, such that there are 
clusters v with the first non vanishing deltas only 
if h v > h; so these clusters give no problems (one 
can even choose the functions \ so that h = 0). Let 
us consider now the case in which the second deltas 
are non vanishing. Note that \(Jhvl~ hv I < lf^7^&l 
< 7 (fc-M(i-O(<0) < 7 (fc-M(V2) as \ ah , \ < jh. h * 

Let us consider vertices v\, V2, ordered so that 
vi > V2 > ... > vi ; then there exists at least a non 
diagonal propagator with scale h m and 

\<J Kil - h ^ I < 7 (*-fc«i)(V2) < 7 (fc„J-fc«/)(V2) 

7 (^-V)(V2) < JK'-KWVJK^-K^im^ 



so that l\ veV2 KJ7-^ < n, e y b 7- (1/2)( ' l ^' v) , if 
Vb are the clusters not verifing the second deltas. 
At the end 

[ J| 7 -[^(P,) + ,(P,)](^-V) ][ ]-J \ ahv \ 7 -^} < 

v VGV2 



(l/2)[D(P v )+z(P v )](h v -h v ,)] 



with D(P V ) + z(P v ) > 0, for any v. Hence a C n e n 
bound follows for |V" fe n |. 
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